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Abstract 

We predict the rate of deep inelastic scattering (DIS) events containing an identified 
forward 7r° that is expected in the experiments at the HERA electron-proton collider, 
in order to see if this process can be used as an indicator of the underlying small x 
dynamics. We determine the background due to deep inelastic events containing forward 
photons which are fragments of the forward jet. We compare the DIS + ir° cross section 
with that of the DIS + parent jet process. 



1. Introduction 



The behaviour of the proton structure function F 2 (x, Q 2 ) at small x reflects the behaviour of the 
gluon distribution, since the gluon is by far the dominant parton in this regime. Perturbative 
QCD does not predict the absolute value of the parton distributions, but rather determines 
how they vary from a given input. If, for example, we are given initial distributions at some 
scale Qq, then the DGLAP p| evolution equations enable us to determine the distributions 
at higher Q 2 . DGLAP evolution resums the leading as ln(Q 2 /Ql) terms. At sufficiently high 
electron-proton centre-of-mass energy ^Js we encounter a second large variable?] 1/x ~ s/Q 2 , 
and we must resum the leading as In 1/x contributions. At leading order the resummation is 
accomplished by the BFKL equation for the (unintegrated) gluon distribution. The solution 
of the equation leads to a singular x~ x small x behaviour of the gluon distribution, where 
A = (3as/7r)41n2 for fixed as and A ~ 0.5 if a reasonable prescription for the running of as 
is assumed and for the treatment of the infrared region 0. The x- x behaviour of the BFKL 
gluon feeds through, via the fc^-factorization theorem Qj, into the small x behaviour of the 
structure function F 2 . Of course in practice we should not expect such a dramatic growth with 
decreasing x, since subleading effects are expected to suppress the effective value of A in the 
HERA domain g]. 

However, it is difficult to identify the presence of the aglnl/rc terms in the measurements 
of F 2 at HERA even though the data do show a steep rise with decreasing x. In fact the rise 
in the latest precise HI and ZEUS measurements || can be well described by next-to-leading 
order DGLAP evolution down to Q 2 ~ 2 GeV 2 and x ~ 10~ 5 . The problem in identifying the 
underlying small x dynamics is due to the parametric freedom that we have in specifying the 
initial parton distributions. For instance for a non-singular gluon input we can increase the 
steepness of F 2 with decreasing x by simply reducing Qq and increasing the DGLAP evolution 
length, \n(Q 2 / Ql) . Alternatively, we could use (as in the global parton analyses [7], §) a singular 
input form xg(x, Ql) ~ x~ x , with A chosen to fit the data. We conclude that it is difficult to 
isolate In 1/x effects from measurements of F 2 . The observable F 2 is too inclusive. Rather, we 
should explore properties of the final state in deep inelastic scattering (DIS). 

The classic way || to probe the small x behaviour of QCD, which avoids the problem of 
assuming input parton distributions, is to study deep inelastic (x, Q 2 ) events which contain 
an identified forward jet (xj, k 2 T ), see Fig. 1(a). According to BFKL dynamics the differential 
structure function for DIS + jet events has the following small {x/xj) behaviour 



d (In 1/2) dk) T ~ C <^i 



4 / 

9 + g {q + q) 



X 



(1) 



where the normalisation coefficient C is given in refs. |10] 11], |12, The parton distributions 



g, q and q are to be evaluated at (xj, k 2 T ). The relevant kinematic region is where 

: To be precise we define x to be the Bjorken variable x = Q 2 /2p.q where Q 2 = — q 2 , and p and q are the 
four momenta of the proton and virtual photon deep inelastic probe respectively. 
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(i) the jet transverse momentum satisfies k 2 T ~ Q 2 so as to neutralize the DGLAP evolution, 
and is sufficiently large so as to suppress diffusion into the infrared region when we solve 
the BFKL equation at decreasing values of x/xj; 

(ii) the jet longitudinal momentum xjp is as large as is experimentally feasible (and x is as 
small as possible) so as to be able to probe the region of small x/xj. 

For these values of Xj the parton distributions entering ([]]) are well known from the global 
parton analyses and so the observation of DIS + jet events offers the opportunity to expose 
BFKL-type small x dynamics free from the ambiguities associated with the choice of the non- 
perturbative parton input. In other words we are studying small x dynamics by deep inelastic 
scattering off a known parton, rather than off the proton. Experimentally, however, the clean 
identification and kinematic measurement of a forward jet proves to be difficult since we require 
it to be as close to the proton remnants as possible, that is Xj as large as possible. Nevertheless 



experimental studies have been attempted and lead to encouraging results pA \. 

Here we use the improved knowledge of the fragmentation functions to propose that the 
forward jet is identified through the measurement of a single energetic decay product. As 
it turns out the 7r° is the hadron which can be identified in the most forward direction in 
the detectors at HERA. We use the BFKL formalism to predict the DIS + forward 7r° rate. 
The rate will, of course, be suppressed in comparison with the DIS + forward jet rate and it 
is an experimental issue to see if the advantages of single particle detection as compared to 
identification of the (parent) jet can compensate for the loss of signal. 

The outline of the paper is as follows. In section 2 we present the QCD formalism required 
to calculate the cross section for the deep inelastic + forward n° process. Then in section 3 we 
discuss the experimental cuts which we impose and give our numerical predictions for the DIS 
+ 7T° cross section. Section 4 contains a discussion. 



2. The DIS + forward 7r° cross section 



First we recall the derivation of the cross section for the deep inelastic + jet process depicted 
in Fig. 1(a), which also shows the variables used. The differential cross section is given by [|K 



dk 2 — 

dk > T xQ* 



(i-y) 



3F 2 



1 2 dF T 

+ ~y 



dxjdk 2 T 2 dxjdk 2 T 



dxdQ 2 

where the differential structure functions have the following leading small x/xj form 



(2) 
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fa (xj, kjj^j $j I , kj T , Q 

n \X<i 



(3) 



2 



for i = T,L. Assuming t-channel pole dominance the sum over the parton distributions is given 
by 

T,f a = 9 + Uq + q)- (4) 

a " 

Recall that these parton distributions are to be evaluated at (xj, k 2 T ) where they are well-known 
from the global analyses, so there are no ambiguities arising from a non-perturbative input. 

The functions <&i(x/xj, k? T , Q 2 ) in (|3|) describe the virtual 7 + virtual gluon fusion process 
including the ladder formed from the gluon chain of Fig. 1(a). They can be obtained by solving 
the BFKL equation 



^i(z,kT,Q 2 



<S>\°\z,k 2 T ,Q 2 ) + 



+ 



3a s 2 

thrp 
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1 dz' dk% 



z Z' JO 



k'S, Q 2 ) - $i(z', k 2 } Q 2 ) ^(z', k 2 , Q 2 



(5) 



The inhomogeneous or driving terms $^ correspond to the sum of the quark box and crossed- 
box contributions. For small z we have 



4°\z,k 2 T ,Q 2 ) » *?>{z = 0,%,C?) = Or). 



(°)^2 



(6) 



We evaluate the by expanding the four momentum in terms of the basic light-like four 
momenta p and q' = q + xp. For example, the quark momentum k in the box (see Fig. 1(a)) 
has the Sudakov decomposition 

k = ap — /3q' + K/p. 

We carry out the integration over the box diagrams, subject to the quark mass-shell constraints, 
and find 



q 



dp / d 2 n T p 2 + (1 - pf 



D\D 2 



^\k 2 T ,Q 2 ) = 2^ el ^ as £ dp J d 2 n T p\l-P) 2 (± - _i 

where the denominators Di are of the form 

D 1 = 4 + p{l - p) Q 2 



Do 



(7) 



(8) 



Do 



[k t - k T ) 2 + P(l - P) Q 2 



assuming massless quarks. 
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If the QCD coupling as is fixed we can solve the BFKL equation (5) and obtain an analytic 
expression for the leading small z behaviour of the solution. Omitting the Gaussian diffusion 
factor in In {k^/Q 2 ) we find 



$ T (z,k 2 T ,Q 2 



$l(z, k 2 T ,Q 2 



512 

2 



-a P +l 



ln(l/z) L 



1 + O 



\n{l/z) 



(9) 
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where the Riemann zeta function ((3) = 1.202 and the BFKL intercept 

12a s 



OCp 



In 2. 



7T 



(10) 



Here, however, we follow the approach of []T2j and allow the coupling as to run. This means 
that we must numerically solve the BFKL equations for 



Hi(z, krp, Q ) — 
We use the differential form of the equations, 



2>a s {k\ 



71 



Hz, k 2 T ,Q 2 ). 



(11) 



d\n{l/z) 



3ag(fcr) t 2 



°° dk$ 



Hi(z, k^-, Q 2 ) H{(z, kji, Q 2 



subject to the boundary conditions 

Hi(z = z ,k^,Q 2 ) 



I V 2 - h 2 I 



H?\k 2 T ,Q 2 ). 



(12) 



(13) 



For the lower limit on the transverse momentum integration we choose k 2 , = 1 GeV 2 . We start 
from the "box" expressions, (0), for H- ' at z — 0.1 and solve (0) to obtain Hi (and $j) for 
z < z . In this way we predict the cross section for DIS + jet production from equations (fj) 
and (g). 

Next let us consider the process where the forward jet fragments into 7r°'s as shown schemat- 
ically in Fig. 1(b). We are looking at the case where the it is collinear with the parent quark 
jet. This means that if the tt° carries a fraction x„ of the proton's longitudinal momentum, 
then it carries a fraction z = x n /xj (0 < z < 1) of the parent jet's longitudinal momentum 
and its transverse momentum k n j< = zk^. In order to calculate the cross section for DIS + 
7T° production we have to convolute the DIS + jet cross section with the it fragmentation 
functions. We obtain 



da* 



dx n dk w T 



dz \ dxi 



dkjrp 



da n 



Df{z,k 2 nT 



2 



da n 



-Df(z,kl T ) 



<i \dxjdk 2 T 
x5 (x-x — zxj) 5 (k n T 



x 



zk jT ) (14) 
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where the sum over q runs over all quark and antiquark flavours. The partonic differential 
cross sections can be obtained from (0) and ([3]) by substituting for the sum over the parton 
distributions Y, a fa either the gluon distribution g or the quark or antiquark distribution |g or 
|g respectively. In analogy to choosing z = 0.1 in flI3|) we impose the constraint xjx-^ < 0.1, 
i.e. x/xj < 0.1 since x n < Xj, on the Xj integration here. The functions Dg°(z,kl T ) and 
Dy° (z, k^ T ) in (|14|) give the probability that a gluon or quark jet fragments into a 7r° carrying a 
fraction z of the parent jet's momentum. We assume that these fragmentation functions satisfy 
leading order DGLAP evolution equations. Note that SU(2) isospin symmetry implies that 

Df(z,kl T ) = l(Df(z,kl T ) + Df(z,kl T )) (15) 



for all partons i — q,g. Therefore ( p!4|) describes the average of the cross sections for n + and 
7r~ production. We use the parametrizations of the leading order charged pion fragmentation 
functions obtained by Binnewies et al. |T5| ; their analysis treated light, s, c, and b quarks 
independently for the first time. 



3. Predictions for the DIS + tt° cross section 



We use (14) to calculate the event rate for deep inelastic scattering in which, in addition to 
the scattered electron, the 7r° is measured in the final state. To ensure that the 7r° is really a 
fragment of the forward jet (and does not come from the quark-antiquark pair which form the 
quark box) we require the 7r° to be emitted in the forward hemisphere in the virtual photon- 
proton centre-of-mass frame. If we express the pion four momentum in terms of Sudakov 
variables 

K = x n p + p n q' + k nT (16) 
then the forward hemisphere requirement is 

> (3 n . (17) 
Since the outgoing pion satisfies the on-mass-shell condition k\ = ~ we have 



x k. 



2 



At = - 7f ■ (18) 

X n 



Then (|17|) gives 



xj > x n > ^xkl T /Q 2 . (19) 

We thus have an implicit lower limit on the Xj integration in (|nj), which is generally stronger 
than the condition Xj > lOx imposed on the solution of the BFKL equation. 

Another problem to be taken into account is that at HERA pions can only be detected if 
they are emitted at a large enough angle to the proton beam. This also ensures that there is 
no contamination from pions produced in the proton remnant. We require 

e*p > 9o- (20) 
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In Fig. 2 we show the relation between the pion kinematic variables for different choices of 
the minimum angle 8 Q defined in the HERA frame. We find that pions with large longitudinal 
momentum fraction x n are only emitted at small angles 6^ p . To reach larger x^ for a given 9^ p 
we can measure pions with larger k 2 T but at a depleted event rate. In the same figure we also 
plot the boundary given by the hemisphere cut, fllTD , for x = 6 x 1CT 4 and Q 2 = 20 GeV 2 , 
which acts as a lower limit on the allowed kinematic region. We will use #0 = 5° for the main 
presentation of our results (although later we compare the predictions with those obtained with 

= 7°). 

Now we are in the position to give numerical predictions for the cross section for the DIS 
+ 7T° production using (|HD and implementing the cuts that we just discussed. Recall that it 
follows from fll~5] ) that the cross section for 7r° production equals the average of the cross sections 
for 7r + and 7r~ production. Therefore the results we will show in the following multiplied by 
a factor of 2 will be valid for charged pion production. Throughout the analysis we assumed 
three flavours of massless quarks. In Fig. 3 we plot the x dependence of this cross section 
integrated over bins of size Ax = 2.10~ 4 and AQ 2 = 10 GeV 2 for three different Q 2 bins, 
namely 20-30, 30-40, 40-50 GeV 2 . Here we required that x n > 0.05 and 3 < k nT < 10 GeV 
and used the fragmentation functions |D| at scale k 2 T . We compare the results obtained when 
BFKL small x resummation is included with the case when gluon radiation is neglected. In 
the first case the strong x dependence of the cross section is driven by the small z behaviour of 
the and therefore there is a strong enhanced increase with decreasing x. For example, if we 
compare the cross section for x ~ 5 x 10 -4 in the two cases, we find that the results are about 
a factor 7 larger when the BKFL resummation is included than when it is neglected. This 
enhancement is the signature of BFKL soft gluon resummation. In fact the BFKL behaviour 
should be identified via the shape in x rather than the value of the cross section, since the 
latter is subject to normalisation uncertainties [|I2" |. In Fig. 4 we show the cross section (in fb), 
for the same cuts as in Fig. 3, in various bins in x and Q 2 which are accessible at HERA. We 
find that the cross section drops off very rapidly with Q 2 which means that we can reach the 
highest values for the bins with 10 < Q 2 < 20 GeV 2 and x very small. 

Of course the DIS + tt° cross section depends on the values chosen for the cuts. In table 

1 we show the effect of changing the limits on the k n T integration. Since the cross section 
decreases with increasing k n T it is more sensitive to the lower limit on the k n T integration than 
to the upper limit. 
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Table 1: The DIS + tt cross section in the bin 1(T 3 < x < 1.2 X 10" 3 , 20 < Q 2 < 30 GeV 2 as 
calculated in Fig. 4, but for different choices of the limits of the integration over the transverse 
momentum of the 7T°. 



KT,min [GeV] 


KT,max [GeV] 


a [pb] 


3 


8 


0.23 


3 


10 


0.26 


5 


10 


0.18 



To obtain the results shown in Figs. 3 and 4 we used the pion fragmentation functions at 
scale k 2 T . In table 2 we show the cross section for the deep inelastic + 7r° process in the bin 
10 -3 < x < 1.2 x 10~ 3 , 20 < Q 2 < 30 GeV 2 calculated imposing the same constraints and 
including BFKL soft gluon resummation but evaluating the pion fragmentation functions at 
the scales \k 2 T) k 2 T and 2k 2 T . The values demonstrate the scale ambiguity in the prediction 
of the cross section. 

Table 2: The DIS + tt° cross section in the bin 10~ 3 < x < 1.2 x 10~ 3 , 20 < Q 2 < 30 GeV 2 
calculated imposing the same cuts as for Fig. 4, but evaluating the fragmentation functions at 
the three different scales |fc 2 T , k\ T and 2fc 2 T . 



fragmentation scale 


a [pb] 


i p 


0.31 


k 2 


0.26 


?k 2 


0.23 



Since 7r°'s are measured through their decay into two photons there is a background from 
events in which the parent jet fragments into a photon which is emitted collinearly, see Fig. 5. 
In analogy to ( |T4] ) the corresponding cross section is given by 



da~, 



dx-fdk^x 



dz 



dkjrp 



D] (z, k 




(21) 



We estimated this background using the fragmentation functions of Owens [0 and found that 
it is 1-2 % of the cross section for tt° production. 

Considering the smallness of the background from photons which are fragments of the 
forward jet, a comment on the errors on the calculation of the cross section for pion production 
is due here. From the numerical point of view there is an error from the Monte-Carlo integration 
used to evaluate (|T4"D which is of the order of 5 %. To our knowledge the errors on the pion 
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fragmentation functions are of the order of a few percent for quarks and 30 - 40 % for gluons 
||17|| . Since the dominant contribution to the cross section for pion production comes from the 
fragmentation of gluons we expect these errors on the fragmentation functions to result in an 
error of at most 25 % on the cross section. The parametrizations of the fragmentation functions 
describe the DGLAP evolution correctly up to 10 % Jy|. We found that our results are more 
sensitive to the normalisation of the fragmentation functions than to their shape. 

Finally let us compare our predictions for DIS + forward tt° production as shown in Fig. 4, 
with the corresponding cross sections for the DIS + forward jet events, the process originally 
proposed by Mueller as the probe of small x dynamics. In order to quantify the suppression 
due to the fragmentation of the jet into the 7r° we integrate the DIS + jet differential structure 
functions given in (|3|) over the same domains of Xj and kj T that we used for x n and k 2 T for 
the DIS + 7T° predictions. To be precise we integrate over the region 3 < kj T < 10 GeV and 
9j P > 5° with a hemisphere cut for the jet in analogy to fllTD , that is Xj > (3j. The upper 
and lower numbers in Fig. 6 compare the DIS + 7T° with the DIS + jet cross section in the 
various bins of x and Q 2 . We see that the fragmentation of the forward jet into a tt° meson 
costs a factor of about 40 in the suppression in event rate. Whether this loss of event rate is 
compensated by the advantage of identifying a forward 7r° as compared to a jet (adjacent to 
the proton remnants) is an experimental question. Table 3 offers a guide to the possible gain 
using the 7r° signal. For instance if we were able to identify 7r° mesons down to 5° in angle and 
5 GeV in kx with the same accuracy as jets down to 7° in angle and 7 GeV in hp then we would 
gain back a factor of 4 0. Moreover if we were to add in the DIS + forward tc^ signal then we 
gain an extra factor of 3. Table 3 also shows that in the HERA regime, where we need to take 
Xj sufficiently large (say Xj > 0.05) to make x/xj small, the low fey events are kinematically 
forbidden by the cuts. For example for 8 Q = 7° we find that kr > 5.0 GeV, while for 6 = 5° 
we have kr > 3.6 GeV. 



4. Conclusion 

In principle, the DIS + jet measurement should be an excellent way of identifying the BFKL 
soft gluon resummation effects at HERA. It turns out, however, that it is experimentally quite 
difficult to measure a forward jet so close to the proton remnants. We therefore suggested 
studying the fragmentation of this forward jet into a single energetic decay product, the n°. 
This should be easier to measure. (The DIS + n° signal can, of course, be supplemented by also 
observing jet fragmentation into ir mesons). We found that when we include BFKL dynamics 
in the calculation of the cross section it leads to the characteristic steep rise with decreasing 

2 Since x^ — zxj with (z) ~ 0.4, our choice of cut-off on x n , that is x^ > 0.05, is probably too conservative |lq |. 
If we were to take x„ > 0.02 then the tt° production rates are all increased. For example if we choose Q 2 ~ 25 
GeV 2 then for x = 5.10~ 4 and 10 -3 the results are enhanced by additional factors of 6 and 5 respectively. Also 
the recent DIS + jet data |l!| lie about a factor of 1.7 above the values that would be obtained from our choice 
of input, so the DIS + tt° rates are expected to be further enhanced by such a factor. 
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Table 3: The DIS + tt° and DIS + jet cross sections in the bin 1(T 3 < x < 1.2 x 1(T 3 , 
20 < Q 2 < 30 GeV 2 as calculated for Fig. 6, but integrated over domains with different choices 
of the minimum angle 9 between the proton and the n° or the jet, and of the minimum 
transverse momentum &/r,mira of the 7r° or the jet. 



k T ,min [GeV] 


#0 


o> [pb] 


Oj [pb] 


3 


5° 


0.26 


10.3 


3.5 


5° 


0.26 


10.3 


5 


5° 


0.18 


8.0 


7 


5° 


0.07 


3.7 


3.5 


7° 


0.08 


3.4 


5 


7° 


0.08 


3.4 


7 


7° 


0.04 


2.0 



x. The disadvantage of using the DIS + tt process is that the event rate is lower than for DIS 
+ forward jet. We quantified the suppression which arises from this jet to ir fragmentation 
(see also footnote 2). It is an experimental question as to whether the loss of event rate can 
be compensated by the more forward domain accessible for n detection and the more accurate 
measurement of the kinematic variables possible for pions as opposed to jets. We presented 
sample results for different acceptance cuts to help provide an answer. 

Since 7r°'s are measured via the two photon decay, there is a background to the deep-inelastic 
+ 7T° measurement from events in which the parent jet fragments into a photon which is being 
emitted collinearly to the jet. We found that this background is about 1-2 %. We conclude 
that deep-inelastic + pion events should be a good way of probing small x dynamics at HERA. 
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Figure Captions 



Fig. 1 Diagrammatic representation of (a) a deep inelastic + forward jet event, and (b) a deep 
inelastic (x, Q 2 ) + identified forward tt° (x n , k n T) event. 

Fig. 2 The relation between the tt° kinematic variables for DIS + 7r° events with x = 6 x 10 -4 and 
Q 2 = 20 GeV 2 for various choices of the angle 9 in ©. In the HERA (27.6 x 820 GeV) 
laboratory frame the pion angle 9 wp to the proton direction is not uniquely specified by 
(x, Q 2 ;x n , k 2 T ). Varying the remaining azimuthal angle transforms the lines of constant 
9 np into narrow bands in the x n , k 2 T plane. Here we averaged over the azimuthal degree 
of freedom. The plot is insensitive to variations of x, Q 2 over their relevant intervals. 
The continuous lines are the upper boundary on the allowed kinematic region for different 
choices of 6 . The dashed line represents the lower boundary given by the hemisphere 
cut, (13), for x = 6 x 10" 4 and Q 2 = 20 GeV 2 . 

Fig. 3 The cross section, (a) in pb, for deep inelastic + tt° events integrated over bins of size 
Ax = 2 x 10~ 4 , AQ 2 = 10 GeV 2 which are accessible at HERA for 7r°'s with transverse 
momentum 3 < k n T < 10 GeV where the constaints x n > 0.05, 9 np > 5°, and the 
hemisphere cut, (|17|) , were imposed. The fragmentation functions were evaluated at scale 
k 2 T . The (a) values are plotted at the central x value in each Ax bin and joined by 
straight lines. The x dependence is plotted for three different AQ 2 bins, namely (20,30), 
(30,40) and (40,50) GeV 2 . The continuous curves show (cr) calculated with $j obtained 
from the BFKL equation. The corresponding (a) values calculated neglecting soft gluon 
resummation and just using the quark box approximation = $j ^ are plotted as dashed 
curves. For clarity a dotted vertical lines joins each pair of curves belonging to the same 
AQ 2 bin. 

Fig. 4 The cross section, (a) in fb, for deep inelastic + tt° events in various (Ax, AQ 2 ) bins which 
are accessible at HERA, and integrated over the region 3 < k^T < 10 GeV, 9 np > 5°, 
x m > 0.05, and subject to the hemisphere cut, ([jj]). The fragmentation functions were 
evaluated at scale k 2 T . The values in brackets are the cross sections obtained when 
using only the quark box approximation = Therefore the difference between the 
two numbers shown in one bin is the enhancement due to BFKL soft gluon resummation. 
Recall that (pl|) implies that the results shown for the DIS + 7r° cross section here equal the 
average of the cross sections for n + and 7r~ production. The curves are the boundaries of 
the acceptance regions at HERA given by 8° < 6 e < 172°, E e > 5 GeV and 0.1 < y < 0.9. 

Fig. 5 Diagrammatic representation of the background to deep-inelastic + tt° events arising from 
photons which are fragments of the forward jet. 

Fig. 6 The upper and lower numbers are respectively the DIS + tt° and DIS + jet cross sections 
(in pb) in various bins of x and Q 2 . For the pion the cuts are those given in Fig. 4, and 
exactly the same cuts are used for the forward jet. 
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